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Abstract 

Using effective field theory methods, we discuss the extraction of the mass and width of the scalar 
mesons /o(980) and ao(980) from the finite- volume spectrum in lattice QCD. In particular, it is argued 
that the nature of these states can be studied by invoking twisted boundary conditions, as well as 
investigating the quark mass dependence of the spectrum. 
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1 Introduction 



The scalar sector of the low-energy QCD is controversial. In particular, in the experimental 
spectrum there are too many candidates for the scalar qq nonet. In the phenomenological 
approaches, alternative solutions to this problem have been suggested. One of the possible 
solutions is to treat some of these mesons as tetraquark states (see, e.g., [1-4]). Another 
suggestion is that a (980) and / (980) are to be considered as KK molecules [5-7]. Further, 
in Refs. [8] these states were described as a combination of a bare pole and the rescattering 
contribution. The investigations carried out within the framework of QCD sum rules are, in 
particular, indicative of the non-gg nature of ao(980) [9]. In the Jiilich meson-exchange model, 
/o(980) appears to be a bound KK state, whereas ao(980) is a dynamically generated threshold 
effect [10]. In the view of the above controversial identifications we wish to stress that all these, 
in general, are model-dependent and can not be unambiguously interpreted in quantum field 
theory. However, in case when the states are very close to some 2-particle threshold (as it is 
indeed the case with ao(980) and /o(980)), it is possible to make a model-independent statement, 
whether these resonances are molecular states or not. The "compositeness criterium," which 
is applied here, was first introduced by Weinberg [11]. This approach is related to the "pole 
counting" method, considered in Refs. [12, 13]. The above methods were used, e.g., in Refs. [14- 
18] to study the nature of ao(980) and / (980) resonances. In particular, in Ref. [16], the 
position of the S'-matrix poles in the vicinity of the KK threshold in the scalar sector of QCD 
is expressed through the so-called Flatte parameters, which describe a resonance located in 
the vicinity of a 2-particle threshold and which are in principle measurable in the scattering 
experiments. Note, however that the compositeness criterium (or the pole counting method) is 
designed to distinguish a near threshold molecular state from a tightly bound system of quarks. 
The question about the precise nature of this quark compound (qq state, or tetraquark, or a 
glueball, or something else), can not be resolved by this criterium. 

It is often stated that the study of the scalar spectrum in lattice QCD can eventually lead 
to the understanding of the nature of these states. Indeed, recent years have seen considerable 
activity, concerning the calculation of the spectrum of the scalar mesons on the lattice (see, 
e.g., [19-25]). However, alone the calculations of the excited spectrum do not answer the 
question. Additional criteria are usually applied. It should be noted that, as compared to 
the phenomenological approaches, lattice QCD has in general more tools at its disposal which 
can be exploited in order to separate the exotic states from the conventional qq spectrum. We 
mention, as one example, the method of hybrid boundary conditions (HBC) [21], which is used 
to distinguish the scattering states from the tightly bound quark-antiquark systems. Another 
example is given by the calculations of the spectrum of the q 2 q 2 mesons, which were done in 
the quenched approximation and in the absence of qq annihilation diagrams [19]. In the latter 
paper it has been argued that, due to the above approximations, and due to the fact that 
the quark masses used in the simulations were rather high, mixing of the qq channel to the 
tetraquark states is suppressed, so that the observed spectrum can be readily attributed to the 
latter. However, even if these arguments might look intuitively plausible, they explicitly refer 
to certain approximations and, for this reason, can not be regarded completely consistent. It 
is evident that, in order to be able to systematically study the nature of the scalar resonances, 
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one should put the existing methods under the renewed scrutiny and look for rigorous criteria, 
which will not be based on the malevolent modifications of the original theory. 

In the present paper we make an attempt to answer the question, how the observables 
of the low-lying scalar mesons a (980) and / (980) can be (at least in principle) determined 
from lattice QCD simulations. It is clear that, due to the proximity of the inelastic threshold, 
finite-volume effects should be very important and could significantly distort the structure of 
the energy levels. Although the finite-volume corrections have been considered in partially 
quenched ChPT at one loop (see, e.g. [25]), a systematic investigation of the problem, to the 
best of our knowledge, is still lacking. The present paper, in particular, intends to fill this gap. 

In addition, we discuss, which conclusions (if any) about the nature of the scalar resonances 
ao(980) and /o(980) can be drawn in a model-independent fashion from the calculations on 
the lattice. Namely, we reformulate the criterium of Refs. [11-14] for energy spectrum of 
lattice QCD in a finite box, whose volume-dependence can be studied within the Liischer 
framework [26]. In order to achieve the goal, using the so-called twisted boundary condition [27] 
has proven to be advantageous. We further investigate the relation of our approach to the HBC 
method of Ref. [21]. 

Further, we discuss a criterium, which can be used to distinguish between the qq mesons 
and the tetraquarks (but which does not distinguish between the tightly bound tetraquarks and 
the KK molecules). The criterium is based on the study of the strangeness content of these 
states. By using Feynman-Hellman theorem, this quantity can be related to the quark mass 
dependence of the exotic state masses and thus can be measured on the lattice. 

The layout of the paper is as follows. In section 2 we briefly review the phenomenological 
determination of the position of a (980) and / (980) poles in the complex plane and discuss the 
"pole counting" method. The generalization of this method to a finite volume is considered in 
section 3. In section 4 we consider the strangeness content of the scalar mesons and formulate 
a criterium for the tetraquark candidates. Section 5 contains our conclusions. 

2 Phenomenological analysis of the tttt scattering ampli- 
tude near KK threshold 

In the following, we restrict ourselves to the discussion of the / (980). The case of the ao(980) 
can be considered along a similar path. 

It is a well-known fact that tttt scattering below KK threshold is almost elastic (the in- 
elasticity parameter is close to unity in this region). For this reason, in the vicinity of the 
KK threshold, where the /o(980) resonance is located, it is convenient to parameterize the tttt 
scattering amplitude in terms of the coupled-channel if -matrix with the following 2-particle 
channels: "l"=7T7r and U 2"=KK. The coupled- channel T-matrix in the S-wave obeys the 
equation 

Tij(s) = Kjj(s) + y]K in (s) iq n {s) T nj (s) , i,j,n = 1,2, (1) 

n 

where qi(s) = ' s/A — M% + iO, (^(s) = a/s/4 — M\ + iO and s denotes the pertinent Man- 
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delstam variable. Note that a similar equation has been used in our treatment of the KN 
scattering length [28]. However, in difference to Refs. [28,29], we do not refer to the non- 
relativistic effective theory in the derivation of Eq. (1). Consequently, the functions Kij(s) 
are no more assumed to be low-energy polynomials in the whole interval between tttt and KK 
thresholds. In principle, Kij(s) has a small imaginary part above 4-7T threshold. However, in 
the following, we shall neglect this effect. 

The S-wave irir scattering amplitude is given by 

Tn(s)= zi . , w 11 = K n + q \ g, . 2 

In the literature, different phenomenological parameterizations of the .fT-matrix are used. We 
distinguish between the parameterizations which have a pre-existing real pole(s) in the vicinity 
of the KK threshold (see, e.g., Ref. [30]) and those which are regular in this region (e.g., [6, 
31]). In the latter case, the i\-matrix elements can be expanded in Taylor series Kij(s) = 
K^' +<z!( s ) Kij ; +0(?2)> whereas in the former, an additional pole term should be also included. 
The location of the S'-matrix pole(s) in either case is uniquely determined by the behavior of 
the i^-matrix in the threshold region. This defines the two-step strategy in the study of scalar 
mesons. In particular, we shall demonstrate below that, measuring the energy spectrum in a 
finite volume, one may uniquely determine the i^-matrix elements that amounts to measuring, 
for instance, the coefficients K^\K^ (this statement is a generalization of Liischer method 
to the multichannel case). At the next stage, continuing the K- matrix into the complex plane 
by using the Taylor (Laurent) expansion, one finds the location of the S'-matrix poles in the 
threshold region from the secular equation 

1 - iq 1 (s)K 11 (s) - iq 2 {s)K 22 {s) - q 1 (s)q 2 (s)(K 11 (s)K 22 (s) - K l2 {sf ) = . (3) 

Assuming further that the quantity w±i in Eq. (2) can be expanded in Taylor series in the 
variable q 2 , one arrives at the so-called Flatte parameterization [32] which modifies the usual 
Breit-Wigner parameterization when a nearby threshold is present 

m const 

T n = -5 1 I 7^ - 5 I I 773 I 7T7~J2Y\ ' ( 4 ) 



m R — s 



0{q 2 ) - im R (g v qi + g K 92(a) + O(gf)) 



where q\ = qi(s) 



Ktf\Kjj. However, it turns out that, for the known phenomenological parameterizations, the 



and the parameters rriR, g„, gK can be expressed through the quantities 



=4M 2 



K 



Taylor expansion of has a very small radius of convergence. For this reason, it is safer to 
use the secular equation, Eq.3 to find the location of the poles. 

The compositeness criterium [11-14] states that, for the molecular states, only one of the 
poles is located in the vicinity of threshold, whereas the nonmolecular state correspond to a 
pair of poles, both lying in the proximity of the threshold. In terms of Flatte parameters, the 
situation gx S> g w corresponds to the KK molecular state and vice versa. 

To summarize, it is clear that the quantities to be measured in the lattice simulations are the 
i^-matrix elements in the vicinity of the KK threshold. These, in turn, determine the position 
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of the poles in the S'-matrix, that is, the energy and the width of /o(980). Consequently one 
can answer the question, whether / (980) is a molecular state or not. If the answer is negative 
this framework does not allow to further elaborate on the structure of this state. For that the 
method discussed in Section 4 has to be used. 

3 i^-matrix formalism in a finite volume 

In Ref. [33] it is shown that the elastic scattering length can be extracted from lattice data, 
studying the volume dependence of the ground state energy. This result is readily obtained 
by Taylor expansion of a more general formula, which enables one to evaluate the elastic 
scattering phase from the finite- volume energy spectrum [26]. At threshold, this scattering 
phase is determined by the effective-range expansion parameters (scattering length, effective 
range, etc), which are thus also obtained from the analysis of the lattice data. Moreover, from 
the phase shift determined on the lattice, one may (in principle) extract the position and width 
of the elastic resonances. A further generalization of the approach allows one to address the 
measurement of the resonance formfactors [34]. For an alternative method to directly extract 
the resonance position in the complex plane from the measured two-point function at finite 
times, see Ref. [35]. 

In the literature there have been attempts to generalize Liischer approach to the case of 
the inelastic scattering (see, e.g. [28,36]). In particular, in Ref. [28], the use of a volume- 
dependent spectrum for the determination of the i^-matrix elements, which are real quantities 
was proposed. The (complex) scattering length can be then expressed through these i^-matrix 
elements. 

The main difference between the approach, Ref. [28] and the present one consists in the 
fact that in the former only the periodic boundary conditions have been used. For this reason, 
it was not possible to determine all three quantities Kij, i,j = 1,2 only from the data taken 
exactly at one energy. Although the effects, coming from the momentum dependence of 
are power-suppressed at large volumes, they still can represent a source of error, if the volume 
is not very large. In the present paper we show that using twisted boundary conditions allows 
one to circumvent this problem. 

Before formulating the approach in a finite volume, few remarks are in order. 

i) The equation, which determines the finite-volume spectrum, has been obtained in Refs. [28, 
29] by using the non-relativistic effective theory. This is the easiest and the most trans- 
parent way of derivation which, however, implicitly assumes that the potential (or the 
.ff -matrix) is a low-energy polynomial. In other words, it is assumed that the Taylor 
expansion of the potential in the momentum space in powers of the relative 3-momenta 
converges for all energies of interest. 

The above requirement is certainly too restrictive, if one considers tttt scattering up to 
1 GeV. On the other hand, this requirement is also superfluous. What is required, is that 
the effective potentials, obtained by the 3-dimensional reduction of the coupled-channel 
mi — KK Bethe-Salpeter equation, are volume-independent up to exponential corrections. 
In fact, this is the path of reasoning, adopted in the original papers by Liischer [26,33]. 



5 



The generalization of the arguments of Refs. [26,33] to the coupled channel scattering 
is relatively straightforward. To this end, first consider the (fictitious) situation where 
2M V > Mk- Then, 2-pion and kaon-antikaon states are two states with the lowest energy. 
The 2-particle irreducible Bethe-Salpeter kernels Uij(P;p,q), i,j = 1,2 in the infinite 
volume are analytic functions of the center-of-mass energy Pq in the range 2M n < Pq < 
4M ff . Then, the regular summation theorem [33] implies that the finite- volume corrections 
to these kernels at large volumes vanish faster than any inverse power of the volume. 

Finally, one should perform analytic continuation in the pion mass. If both pion and kaon 
masses are taken physical, inelastic thresholds move below the KK threshold. Strictly 
speaking, it is not true any more that corrections to the Bethe-Salpeter kernels are ex- 
ponentially suppressed in the energy interval we are considering. However, as already 
mentioned, the coupling to the inelastic channels is extremely weak: the inelasticity pa- 
rameter 7] > 0.98 below KK threshold. This means that the coupling to the inelastic 
channels can significantly influence the spectrum obtained in the absence of these chan- 
nels, if the two-particle and multiparticle energy levels in a given volume accidentally come 
very close to each other. Neglecting this possibility, we expect that inelastic channels ef- 
fectively decouple and the kernels can be considered to be almost volume-independent. 
We shall use this assumption below. 

ii) On the cubic lattice the rotational symmetry is broken down and mixing of all partial 
waves occur. For the problem in question this mixing, however, is strongly suppressed. In 
order to see this, let us recall that if the cubic symmetry is not broken, the S-waves mix 
with the partial waves with the orbital momentum I > 4 [26]. Below 1 GeV, the partial 
wave with I = 4 is however very small (see, e.g. [37]). Using Eq. (6.15) of Ref. [26] and 
the parameterization of the G-waves from Ref. [37], we have estimated the correction term 
arising from the mixing. In the region of interest it does not exceed a few percent. In the 
following, we neglect the mixing altogether. 

To summarize, with the above assumptions, the generalization of Eq. (1) to the finite- volume 
case is straightforward, and the result is similar in form to that given in Ref. [28]: 

Tij{s) = Kij{s) + — — ^K in {s) Z 00 {l,kl)T nj (s) , i,j,n= 1,2, (5) 

* n 

q n (s)L 

where k n (s) = — and Zoo(l;& 2 ) stands for Liischer zeta-function [26,33]. 

2ir 

The equation (5) implies periodic boundary conditions on both pion and kaon fields. Below, 
we explore the possibility of using the so-called twisted boundary conditions, which we opt 
to impose only on the strange quark field, whereas u- and ci-quarks obey periodic boundary 
conditions 

u(x + Le^ = w(x) , d(x + Le^ = d(x) , s(x + Lei) = e ie s(x) , < 9 < 2ir . (6) 

In this relation, e«, % = 1,2,3 denote unit vectors along the lattice axes, and V = L 3 is the 
spatial volume of the lattice. The vacuum angle theta is chosen the same in all spatial directions, 
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L (in units of M 



Figure 1: Energy levels E n , n = 1,2, ■•• of the two-pion states in a finite box. Solid lines 
correspond to the periodic boundary conditions (9 = 0), and the dashed lines to the antiperiodic 
boundary conditions for the s-quark (9 = it). Horizontal dashed line depicts the KK threshold. 
Different parameterizations of the i^-matrix from Refs. [6,30,31] have been used. 



in order to avoid the breaking of the cubic symmetry. This choice, however, may be changed, 
if needed. 

In the effective theory, the angle theta will appear in the boundary condition for the kaon 
field and not for the pion field. 

i^(x + Lei) = e Tie K ± (x) , K°(x + Le { ) = e~ i6 K°(x) , K°(x + Le t ) = e i0 K°(x) . (7) 

/ Q Q Q\ 

If 9 7^ 0, the KK pair at rest has the minimum relative momentum ± ( — , — , — j and the energy 



W 

2\\M\ + —j. In other words, the KK threshold can be moved by adjusting 9. Exactly this 

property makes the twisted boundary conditions particularly useful to study scalar mesons, 
which are located very close to this threshold. 
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Parameterization: Oiler and Oset 





Figure 2: The tangent of the pseudophase in the vicinity of the KK threshold, in case of the 
periodic boundary conditions (6 = 0). The pseudophase was extracted from the excited energy 
levels E n with n = 3,4,5. Different parameterizations have been used. 



Twisted boundary conditions can be straightforwardly implemented in the two-channel 
Liischer equation (5). To this end, it suffices to replace Z 00 (l; k%) — > Zq (1; k%), where 



Z°(l, k\) = lim V — 1 —— , (8) 

whereas Z 00 (l; k\) stays put. 
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The energy levels are given by the secular equation, which is obtained from Eq. (5) 



1 " -^r £oo(l; k\) K n (s) - -L- Z e m {\- k 2 2 ) K 22 (s) 

Z 00 (l; hi) -I- Z 9 (l; A; 2 ) (ifii(s)if 22 (s) - K l2 {sf) = . (9) 



1 v^L~ uu ^"^VEL 

Our aim is to describe a procedure, which enables one to extract the if -matrix elements ifjj 
in the vicinity of threshold from the lattice. We shall illustrate this procedure on the example 
of the synthetic data, which were produced, using different phenomenological parameterizations 
for the inr S-wave phase shift and the Eq. (9) to determine the spectrum for various values of 
9. We have tested the parameterizations, given in Refs. [6,30,31] (these are shortly described 
in appendix A), and produced the spectrum which is shown in Fig. 1. In this figure, the 
spectra at 9 = (solid lines) and 9 = n (dashed lines) are displayed. These spectra show quite 
similar behavior despite the fact that the pertinent if -matrices have very different properties. 
For example, the if -matrix from Ref. [30] has a real pole very close to the KK threshold, 
whereas the if-matrix from Refs. [6, 31] is regular near threshold. As expected, the spectrum 
is almost independent on the twist parameter 9 away from threshold. Maximal variation is 
introduced in the vicinity of the threshold where the rearrangement of the levels occurs: the 
levels with 9 = ir for small values of L are "pushed up" one level high as compared to the case 
with 9 = 0. Consequently, the measurements for different values of 9 provide an independent 
piece of information in the threshold region, where the ^-dependence is maximal. Since we are 
looking for the resonances exactly in this region, twisted boundary conditions can be used to 
fix the location of these resonances. On the other hand, it is also clear from Fig. 1 that the 
attempts to identify separate levels with either the resonance or scattering states are not very 
informative. In particular, tuning the parameters L, 9, one may easily move a single energy 
level above or below the threshold. 

In order to facilitate the extraction of the if-matrix elements from the data, according to 
Ref. [28], it is convenient to define the pseudophase 

tan 5° (ft) = -tan0(fci) , k x = — , tan (j){k x ) = - . (10) 

2vr Z 00 (l;k() 

The {9- dependent) energy spectrum ft = qi(L) can be (in principle) measured on the lattice. 
Consequently, the right-hand side of Eq. (10) (and, hence, the pseudophase) are measurables. 
The physical meaning of the pseudophase is the following: apply Luscher formula to the mea- 
sured energy spectrum, assuming that the scattering is elastic (in our case, this means that only 
7T7T threshold is taken into account). Thus, below the inelastic KK threshold, the pseudophase 
coincides with the usual scattering phase. This is not the case above the inelastic threshold. 
From Eq. (9) one may express the pseudophase through the if-matrix elements as follows: 



/ ^ Z e m (V,kl) K 12 (sf \ 
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Suppose, for a moment, that 9 = and we tuned the box size so that the energy is exactly 
equal to 2Mk (this corresponds to L = L\ in Fig. 1) and measured the pseudophase for this 

box size. Recalling that Z 00 (l; fc|) = 7== 2.514488997 + 0(H) as k\ -> 0, from Eq. (11) 

V47r/;;o 

we readily get 



2 



lim tan5 e (g 



e=o V K\ 2 ' 



;(o) 



This equation gives one relation between three quantities if^. 

We repeat this procedure for a different value L and adjust the parameter 9 so that energy 
of the measured level is exactly 2Mk again. For example in Fig. 1 this corresponds to L = L 2 
at 9 = 71. After performing three measurements at threshold energy and different values of L, 9, 
we get enough equations to determine all K^' separately. Moreover, there is nothing special 
about the threshold energy: the same procedure can be repeated at ^Js 7^ 2Mk, scanning 
the matrix elements Kij(s) in the vicinity of threshold. In this way, one may e.g., answer the 
question, whether the .fT-matrix contains the pre-existing poles in the threshold region. 

From Figs. 2 and 3 one may conclude that the difference between different parameteriza- 
tions of the i^-matrix is more clearly visible in the pseudophase that in the structure of the 
energy levels which all show a similar behavior. Moreover, it is seen that the behavior of the 
pseudophase changes dramatically when 9 changes from to ir. On the basis of this observation 
one may expect that the equations that relate the matrix elements Kij(s) with the measured 
pseudophases at the same energy are not degenerate and will enable one to neatly extract these 
matrix elements from the data. 

Last but not least, we wish to comment on the relation of the approach suggested in the 
present paper with the method of HBC [21]. We shall do this, adapting the argumentation 
of Refs. [21] to a choice of the boundary conditions used in the present article. Namely, as 
already mentioned above, varying the parameter 9 from to it amounts to floating the KK 
threshold, whereas the tttt threshold stays put. Consequently, one interprets a given state as 
a KK scattering state if it is dragged along by the KK threshold, and as a "genuine quark 
state," if it does not move. On the other hand, from the expression of the pseudophase Eq. (11) 
one may conclude that the energy level in the vicinity of the threshold does not move if the 
matrix element Kyi that describes the coupling of tttt and KK channels, is small. In terms of 
Flatte parameters, this corresponds to a small value of gx, see Eq.4. In this case, a pair of poles 
appears in the S'-matrix near the KK threshold. To summarize, it is seen that the approach 
proposed in the present paper is a natural generalization of the HBC method and enables one 
to extract and to quantify the information about the nature of the resonances in the threshold 
region. 

Finally, we wish to mention that we performed the calculation of the energy levels and 
the pseudophases (for the different values of the twisting angle 9) for the case of the a (980) 
resonance as well, using the parameterization of Ref. [6]. The results turn out to be qualitatively 
similar to the case of the / (980). 
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Figure 3: The same as in Fig. 2, but with antiperiodic boundary conditions 6 = tt. 

4 Strangeness content of the exotic states 

As mentioned in the introduction, different methods are used at present to distinguish tetra- 
quarks from ordinary qq mesons in lattice QCD. Usually, a state is said to be a tetraquark if it 
is seen in the two-point function of the operators qqqq and not seen in the two-point function 
of the operators qq. More precisely, this means that the overlap of the one particle state with 
the state produced from the vacuum by the operator qqqq is much larger than with the state 
produced by the operator qq. Furthermore, one may invoke arguments based on the behavior 
of the spectrum in quenched approximation [19]. However, one should take these arguments 
with a grain of salt. For example, the the matrix elements used in the above argumentation 
are scale-dependent and, in addition, depend on the way the operators are constructed from 
the quark fields. It is clear that a mathematically consistent and model-independent criterium 
would be highly desirable. 
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B 


qq 


q 2 q 2 


1 = 0, nonstrange 
1 = 0, strange 

I = 1 


uu+dd „, n 

ss, y B = oo 
us, ds + conj., ub = 2 
ud, du, y B = 


[ud}[ud], yB = 

{su][su]+[sd][sd} _ 
V2 2/B — ^ 

[sdlfwd] + conj., Vb = \ 



Table 1: The wave functions and the strangeness content of gg and g 2 g 2 mesons in the non- 
relativistic quark model with ideal mixing. 



It should be pointed out that, generally speaking, the question whether the composition 
of a particular state is qq or q 2 q 2 , makes sense e.g., in the non-relativistic quark models but 
not in the field theory. In the latter, any operator with appropriate quantum numbers can 
be used as an interpolating field for a given particle. Thus, we have to look for a criterium 
formulated in terms of observable quantities, which in the non-relativistic limit is coherent with 
our understanding of ordinary qq and tetraquark states. We shall use this strategy in the 
following. 

Note first that the quark model wave functions of the states in the non-relativistic limit are 
eigenf unctions of the operator 

i = u,d,s 

with the eigenvalues Ni + Nj, where iVj and Nj denote the number of quarks and antiquarks of 
a flavor % which are contained in this state. 

Consider now the SU(3) nonets with maximal mixing for the qq and q 2 q 2 mesons (the wave 
functions of the q 2 q 2 states are given, e.g., in Ref. [38]). Further, define the strangeness content 
of a state B that belongs to these nonets in a standard manner 

2(B\ss\B) . 
\B\uu + dd\B) 

In the non-relativistic case, in order to calculate ys, one has just to count the number of the 
quarks (antiquarks) of a given species in a state B. Using the wave functions from table 1, we 
may easily evaluate the values of ub- The results are given in the same table. 

As one sees from table 1, the patterns followed by ys are very different for the qq mesons 
and tetraquarks. Note also that, in case of arbitrary mixing, these two nonets can be still 
distinguished due to the fact that ys for / ^ does not depend on the mixing angle. 

Next, we mention that the quantity yB is a well-defined quantity in QCD (it is scale- 
independent) and can be directly evaluated on the lattice e.g., by studying the quark mass 
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dependence of the scalar meson masses and applying Feynman-Hellmann theorem 



Vb = 2 




dM, 



dm 



B 



) 



(15) 



where we have assumed that isospin is conserved: m u = m& = m. 

Now we are in a position to formulate our proposal. The quantity ?/b is a well-defined 
quantity in QCD and can be measured on the lattice. On the other hand, in the non-relativistic 
quark models this quantity clearly distinguishes between the qq and q 2 q 2 states. Therefore we 
can use the measured pattern of the quantity ys for the nonet states in order to define, what we 
understand under tetraquark states within relativistic quantum field theory. In contrast to the 
criteria used in the literature so far, this definition, e.g., does not operate with the quantities 
that are scale-dependent (like the matrix elements of multiquark operators). Note that we took 
advantage here of being able to vary quark masses on the lattice freely. Such a possibility is not 
available in the phenomenological approaches based on the experimental input and the above 
definition will be harder to use there. 

Finally, note that the scalar mesons under consideration are resonances, not stable particles. 
Since the width of these resonances is very small, it is natural to continue using Eqs. (14) and 
(15). For example, Mb in Eq. (15) is to be now understood as the resonance pole position 1 . 
This means that the measured energy levels should be first "purified" with respect to the finite 
volume effects, as described in the previous sections. The method should be applied at the 
final stage, to the resonance poles extracted from the spectrum. If this is not done, the KK 
threshold, which will be moving if quark masses are varied, could strongly influence nearby 
energy levels, and this may result in wrong conclusions about the quark mass dependence of 
the true resonance energies. 

Interestingly, the kaon mass dependence of the scalar mesons can also be used as a signature 
of the molecular picture. For a such a scenario, the leading Fock component of the scalar meson 
wave function has two quarks and two anti-quarks, much like the just discussed tetraquark 
states. However, while these are expected to be compact, a molecule is loosely bound and thus 
spatially extended. The molecular nature leads to a very peculiar kaon mass dependence of the 
molecule as shown in Ref. [39]. The mass of a KK- molecule can be written as M mo \ = Mk + 
M K — e, with e the small binding energy, e <C Mk- Consequently, the kaon mass dependence of 
M mol is expected to be linear with a slope of two, with correction of the order C(r 2 eM mo i), where 
r denotes the range of forces. This is dramatically different from a tetraquark, where the kaon 
mass dependence is generated either from the strange valence or the strange sea quarks and is 
thus expected to depend on m s linearly or as m Z J 2 , generating a leading kaon mass dependence 
as M K or M K , respectively. For a more detailed discussion on this issue, we refer to Ref. [39]. 



The procedure of the extraction of the resonance matrix elements on the lattice is discussed, e.g., in Ref. [34] . 
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5 Conclusions 



The main results of this investigation can be summarized as follows: 

i) Obviously, the measured excited spectrum can not be directly identified with the exper- 
imentally observed scalar mesons. This simple fact becomes crystal clear by looking at 
Fig. 1, where one has freely used the parameters L and 9 to move an energy level below 
or above 2Mk- What defines the energy and the width of a resonance is the position 
of the S'-matrix pole in the complex plane. This position can be determined, extracting 
i^-matrix elements from the spectrum measured on the lattice. 

ii) The procedure of determining i^-matrix elements at the inelastic threshold, which was 
described in the present paper, is the generalization of Liischer's method to the elastic 
scattering length. It is also an improvement of the method described in Ref . [28] . Namely, 
using twisted boundary conditions, as proposed in the present paper, enables one to de- 
termine all three quantities Kij(s) in the vicinity of the inelastic threshold. 

iii) Liischer's approach implies the study of the response of the energy spectrum on the varia- 
tion of the box size L. We have shown that, for certain quantities, studying the dependence 
on the twisting angle 9 may partly substitute studying the volume dependence. In other 
words, one may use the twisting parameter 9 to scan the energy region in the vicinity of 
the KK threshold. 

iv) The study of the L- and ^-dependence of the energy spectrum of scalar mesons as proposed 
in the present paper is beyond any doubt a very demanding enterprise. The authors bear 
no illusion that the whole program can be realized in lattice calculations at physical quark 
masses anytime soon, especially for the /o(980) meson (the situation with ao(980) could 
be slightly better). However, we still find it important to formulate a rigorous way to treat 
the problem in question, which can be used one day. 

v) We show that the use of the twisted boundary conditions allows one to distinguish between 
the loosely bound molecular states and the compact quark compounds. We in addition 
argue that if the latter possibility is realized, the measured values of the strangeness 
content for the different members of the SU (3) nonet allow one to interpret these states 
either as conventional qq states or q 2 q 2 tetraquark states. Thus the measurement of the 
strangeness content on the lattice, which can be achieved by studying the quark mass 
dependence of the resonance energies, enables one to gain detailed information about the 
structure of the scalar mesons. We have also pointed out that the molecular picture can 
be further tested from the measurement of the kaon mass dependence of the mass of the 
scalar mesons - for a molecular state this would be linear with slope two (modulo small 
corrections) . 
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A Different parameterizations of the two-channel 
K- matrix 

A particular parameterization of the A'- matrix elements from the paper Protopopescu et at., 
Ref. [31], which is used in the present paper, is given by 



A 



(m- 1 ; 



Mij(s) 



M. 



(o) 



s t )M, 



(i) 



ZJ 1 



s t = 4M| 



(a.i; 



where the coefficients Ml- take the following values 



'j 





M$ 






n = 


3.38 M w 


2.40 M n 


0.071 M n 


n = 1 








-0.0038 AAr 1 



The scattering matrix in the paper by Oiler and Oset, Ref. [6], is given by a solution of the 
2-channel Bethe-Salpeter equation 



Tij(s) = Vij(s) + ^ V in {s)G n {s)T nj {s) , i, j, n = 1, 2 



(A.2) 



where Vij(s) are the tree- level meson-meson scattering amplitudes, calculated in Chiral Pertur- 
bation Theory 



2F 2 ' ' " ' ^ "" 4A 2 ' ' 22 " 4F 2 ' 

7T 7T 7T 

with F n ~ 93 MeV and the loop functions Gfe(s) are given by 



-M 



3s 



G„ = A/" -1 



1 ^ 



i 



On ln <Tn< a 7? r 



16vT 2 o n w™*l(f 



M=-(%-kJs)- 1 (A.3) 



G* + zg n ( S ). (A.4) 



8tt 2 ~ M n 

In the above expression, a n = yl 1 — + iO, g max = a/A 2 — M\ is the cutoff momentum in 

the loops and w™ ax = a/ M 2 + (g max ) 2 . The cutoff parameter A is chosen to be A ~ 1020 MeV. 
The A-matrix elements are given by a solution of the equation 



(A.5) 
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In analogy to the parameterization by Protopopescu et al, the above i^-matrix elements are 
regular in the vicinity of the KK threshold. The resonance emerges due to the rescattering 
effect. 

In difference to this, the parameterization of the K- matrix in the paper by Au et al.,, Ref. [30] 
contains the pre-existing pole in the vicinity of the KK threshold. This parameterization looks 
as follows 

^W k-^l +g^- 1 )} (A ' 6) 

where s = —0.0162, si = 0.9383, f\ = —0.1659, f 2 = 0.5852 and the values of the coefficients 
c™ are given by 





„n 
c ll 


„n 
c 12 


„n 


n = 


0.4247 


-3.1401 


-2.8447 


n = 1 


-0.5822 


-0.1359 


6.9164 


n = 2 


2.5478 


1.0286 


5.2846 


n = 3 


-1.7387 


-2.3029 


-0.9646 


n = 4 


0.8308 


0.1944 






All dimensionful quantities are given in powers of GeV. 
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